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Abstract 

We present the exact solution of a two-dimensional directed walk model of a drop, or 
half vesicle, confined between two walls, and attached to one wall. This model is also a 
generalisation of a polymer model of steric stabilisation recently investigated. We explore 
the competition between a sticky potential on the two walls and the effect of a pressure-like 
term in the system. We show that a negative pressure ensures the drop/polymer is unaffected 
by confinement when the walls are a macroscopic distance apart. 

1 Introduction 

The study of the behaviour of the boundary between two phases has a long history [HE]- In 
two-dimensions where the boundary is one-dimensional much work has been done to understand 
such behaviour [3]. One type of model that has proved useful involves directed walks in a half 
plane with a Boltzmann weight associated with the area under the walk [JJ. Related to this is 
the study of an enclosed boundary that can be used to model biological membranes [5]. The 
study of lattice vesicles is motivated by biological membranes that consist of lipid bi-layers and 
form shapes that depend on acidity, osmotic pressure and temperature. The basic model in two 
dimensions consists of some type of self-avoiding polygon on a lattice which is weighted according 
to its area and perimeter. The weighting of the area is analogous to an osmotic pressure. Various 
exactly solved cases including directed walk model have been considered [6j El [8] . 

In a parallel development the study of the behaviour of a long linear polymer molecule in 
dilute solution confined between two parallel plates [9] has recently gained momentum with the 
exact solution of various directed walk models |10|. II lj . The phenomena being modelled here are 
the steric stabilisation and sensitised flocculation of colloidal dispersions. The progress made 
by the exact solution is the finding that if one considers a polymer confined between two sticky 
walls, the thermodynamic limit where the polymer is much longer than the distance between 
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the walls is shown to have a different phase structure to the case where is there is one wall only 

ma. 

In this paper we study a directed version of the self-avoiding vesicle model that is a general- 
isation of both the vesicle models studied previously [8] and the recently studied polymer model 
[TO] , It is effectively a half- vesicle, or drop, if considered as a phase boundary. We consider 
directed self-avoiding walks on the square lattice, confined between two lines (y = and y = w) 
whose ends are both attached to one of the lines so that it effectively forms a polygon or rather 
a loop. We add weights for visits of the walk to both the top wall and the bottom wall. We also 
add an (osmotic) pressure-like term that weights the area contained under the loop. 

Here we solve this model to show that introduction of such a negative (osmotic) pressure 
ensures that the two wall scenario behaves in the same way as the one wall case, in contrast to 
the model without any osmotic pressure. 

2 The model 

The directed walk model that we consider is closely related to Dyck paths. Dyck paths are 
directed walks on Z 2 starting at (0, 0) and ending on the line y = 0, which have no vertices with 
negative y-coordinates, and which have steps in the (1,1) and (1,-1) directions. We impose the 
additional geometrical constraint that the paths lie in the slit of width w defined by the lines 
y = and y = w. We refer to Dyck paths that satisfy this slit constraint as loops (see figured]). 



Figure 1: An example of a directed path which is a loop: both ends of the walk are fixed to be on 
the bottom wall. 

Let C w be the set of loops in the slit of width w. We define the generating function of these 
paths as follows: 

L w (z,a,b,p)= z n Wo tt Wt"Wp mW ; (2-1) 

where n(ir) , u(tt) , v(n) and m(n) are the number of edges in the path ir, the number of vertices 
in the line y = (excluding the zeroth vertex), the number of vertices in the line y = w, and 
the area under the walk counted as the number of half-faces of the lattice between the walk and 
the line y = 0, respectively. 

Let £™ be the sets of loops of fixed length n in the slit of width w. The partition function 
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of loops is defined as 

Z n (w;a,b,q) = a «W6«W g ™(*) . ( 2 .2) 

Hence the generating function is related to the partition functions in the standard way, and for 
loops we have 

oo 

L w (z,a,b,p) = '^2z n Z n (w;a,b,p). (2.3) 

n=0 

We define the reduced free energy k(w; a, b, q) for loops for fixed finite w as 

k(w; a, b,p) = lim n~ l log Z n (w; a, b,p) . (2-4) 

n— »oo 

Consider the singularities z c (w; a,b,p) of the generating functions L w (z, a,b,p) closest to 
the origin and on the positive real axis, known as the critical points. Given that the radii of 
convergence of the generating functions are finite, which we shall demonstrate, and since the 
partition functions are positive, by Pringsheim's Theorem (see e.g. Theorem IV. 6 in [12] ) the 
critical points exist and are equal in value to the radii of convergence. Hence the free energies 
exist and one can relate the critical points to the reduced free energies as 

K(w;a,b,p) = -logz c (w;a,b,p). (2.5) 

3 The half-plane 

3.1 Preliminaries 

We first consider the case in which w — > oo, which reduces the problem to the adsorption of 
paths to a wall in the half-plane with osmotic pressure. 

We begin by noting that once w > n for any finite length walk there can no longer be any 
visits to the top surface so 

Cl +1 = Cl for w > n. (3.1) 

Hence we define the sets C 1 ^ as 

£-hp = • (3-2) 

The limit w — > oo can therefore be taken explicitly. 

Also, as a consequence of the above, for all w > n we have v(ir) = for any ir £ £™ . Hence 
the partition function of loops in the half-plane can be defined as 

Z h n p (a,p)= £ a"Wp m W . (3.3) 
We define the generating function of loops in the half plane via the partition function as 

oo 

L h P{z,a,p) = Y J Z n Z^{a,p). (3.4) 

n=0 
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In an analogous way to the slit we define the reduced free energy in the half-plane for loops as 
K hp (a,p) = lim n-HogZ^i^p) = - logz£ p (a,p). (3.5) 

71— >00 

We note that in defining these free energies for the half-plane the thermodynamic limit n — > oo 
is taken after the limit w — > oo: we shall return to this order of limits later. 

3.2 Exact solution for the generating function in the half-plane 

We use a standard decomposition argument to derive a functional equation for the generating 
function L hp (z,a,p) as follows. Except for the zero-step Dyck path with weight 1, every Dyck 
path can be decomposed uniquely into a Dyck path bracketed by a pair of up and down steps, 
followed by another Dyck path. The associated generating functions are apz 2 L hp {pz, l,p) and 
L hp (z, a,p), respectively. This decomposition leads to the functional equation 

L hp (z,a,p) = l + apz 2 L hp ( y pz,l,p)L hp (z,a,p) , (3.6) 

which one can rewrite as 

LhP(2 ' a ' P) = l-apz 2 L^(pz,l,p) ■ (3 - 7) 
By iterating this equation one can find a continued fraction expansion for L hp (z,a,p) as 

L hp (z,a,p) = — . (3.8) 

apz 



p 3 z 2 



^ p 5 z 2 



p 7 z 2 



I- ... 

One can find this solution in the literature |14j . To find a series solution one can use a lin- 
earisation Ansatz, standard for a g-deformed algebraic equation such as Equation (|3.6|) . We 
substitute 

**(M,,) = ^ (3.9, 
into Equation (|3.6|) with a = 1, and find that H(z, q) must satisfy the linear (/-functional equation 

H(pz,p) - H(z,p) -pz 2 H(p 2 z,p) = . (3.10) 

We then solve this linear functional equation using a series in z, 

oo 

H(z,p) = Y^c n (p)z n . (3.11) 

n=0 

This leads to the simple two-term recurrence 

(l-p n )c n -p 2n - 3 c n . 2 = . (3.12) 
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By iteration and using initial conditions cq = 1 and c\ = one finds 

00 2m(m-l) 

ff (*>p) = E ^ 2 2 , (-p* 2 r • (3.13) 

Here (t;q) n = Yl^Z (l — tq k ) is the standard (/-product. Using basic hyper geometric series 
notation [13], we identify H(z,p) = o4>i(— ;0;p 2 , —pz 2 ), where q0i(— ;b;q,t) is given by 



^ (b;q) n {q;q)n 



Via (|3.7|) and (|3.9p we arrive at the final result 



L h P(z,a,p) = H{PZ ' P) 

H(pz,p) — apz 2 H(p 2 z,p) 



aH(z,p) + (1 - a)H(pz,p) 
where H(z,p) = o0i( — ; 0;p 2 , —pz 2 ). 



H{PZ > P) (3.15) 



4 Exact solution for the generating functions in finite width 
strips 

The solution of the strip problem begins by using an argument that builds up configurations 
(uniquely) in a strip of width w + 1 from configurations in a strip of width w. In this way 
a recurrence-functional equation is constructed. Consider configurations of loops in a strip of 
width w (see figure [2]), and focus on the vertices touching the top wall: call these top vertices. 
These vertices contribute a factor b to the Boltzmann weight of the configuration. Now consider 
a zig-zag path (see figure [2]), which is defined as a path of any even length, or one of length 
zero, in a strip of width 1. The generating function of zig-zag paths is 1/(1 — bpz 2 ). Replace 
each of the top vertices in the configuration by any zig-zag path. Since one could choose a single 
vertex as the zig-zag path all configurations that fit in a strip of width w are reproduced. Also, 
the addition of any non-zero length path at any top vertex will result in a new configuration 
of width w + 1 and no more. The inverse process is also well defined and so we can write 
recurrence- functional equations for each of the generating functions. The generating function 
L w (z, a,b,p) satisfies the following functional recurrence: 

Li{z,a,b,p) = ^7—0; (4.1) 

1 — abpz z 

L w (z,a,b,p) = L w _i fz,a, ^ _ b ^ 2 w-i z 2 >p) ■ ( 4 - 2 ) 

We note that the zeroth vertex of the path is weighted 1 and that L w (z, a, b,p) counts the walk 
consisting of a single vertex. We see that 

1 

L 2 (z,a,b,p) = — (4.3) 



apz 2 



1 — bp A z 



3^2 
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Figure 2: TTte construction of configurations in strip width w + 1 /rom configurations in strip 
width w is illustrated. Every vertex touching the upper wall in the strip of width w can be replaced 
by a zig-zag path as shown. 



and 



L 3 (z,a,b,p) = . (4.4) 

apz 

p 3 z 2 

1 



1 — bp 5 z 2 

Hence, by substitution, this functional recurrence generates a finite continued fraction expansion 
of the generating function 

1 

L w (z,a,b,p) = . (4.5) 

apz 



1 - 



p 3 z 2 



p 5 z 2 



p 7 z 2 



„2«)— 3 r 2 

1 - bp 2w - l z 2 

This finite continued fraction can be compared directly to (|3.8p . which is formally obtained by 
taking the limit of w — > oo. Conversely, L w (z,a,b,p) is a finitary version of L hp (z,a,p). 
It is clear that the generating function can also be written as a rational function 

t ( h \ Pw(z,a,b,p) 

L w (z,a,b,p) = —— — — (4.6) 

Q w (z,a,b,p) 

though it does not simply follow to write expressions for these. It does however follow from 
the theories of continued fractions and orthogonal polynomials (see pages 256-257 of Andrews, 
Askey and Roy [15] ) that both the numerator P w and denominator Q w of the generating function 
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satisfy recursions 



P w (z,a,b,p) = < 



1 , w = 1 

1 — bp 3 z 2 , w = 2 

P m _i(^,a,l,p) - bp 2w - 1 z 2 P w ^ 2 (z,a,l,p) w>3, 



and 



Q w (z,a,b,p) = < 



1 — a6pz 2 , w = 1 

1 — apz 2 — bp^z 2 , w = 2 

Q w -i(z,a, l,p) - bp 2w ~ 1 z 2 Q w ^ 2 {z,a, l,p) , w > 3 . 
One can immediately note that 



so that 



P w (z,a,b,p) = Q w (z,Q,b,p) 

T , , x Q w {z,0,b,p) 
L w (z,a,b,p) - 



Q w (z,a,b,p) 

We now form the width generating function for the denominator as 

oo 

W(t,z,a,b,p) = ^] Q w (z,a, b,p)t w , 

w=0 

and find a functional equation for W(t) from the recurrence (|4.8[) as 



W(t,z,a,l,p) - 

oo 

E 

n=0 



(i;p 2 )n 



+ 1 



(4.7) 



(4.8) 



(4.9) 
(4.10) 

(4.11) 



W(t,z,a,b,p) = 

t(l - i)(l - a&pz 2 ) + i 2 (l - apz 2 - op 3 z 2 ) + tW(t, z, a, l,p) - 6p 3 z 2 t 2 W(p 2 t, z, a, l,p) . (4.12) 
One can first solve for W(t, z, a, l,p) by iteration to give 



(4.13) 



It is then straightforward to provide an expression for W(t, z,a,b,p) by substituting into the 
functional equation ()4.12p . To simplify the expressions further, let us rewrite our expression for 
W(t, z, a, l,p) in terms of 

n n(n— l)^,n 

(4.14) 



~ qn(n-l) x n 
<l>(t,q,x) = 2^ 



n=0 



(t,q)r 



Using basic hyper geometric series notation [T3], we identify 4>(t,q,x) = 1^2(0; 0, i; g, x). With 
help of <f)(t, q, x) we can write 



W(t,z,a,l,p) 



1 — apz z 
pz 2 t 



1 



1 — apz 2 



pz 2 t 



0(t,p z , -pzH 2 ) + <f>(t,p z , -p 6 z 2 t 2 ) , (4.15) 
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which leads to 



W(t,z,a,b,p) 



b + t — bt — apz 2 t — bpz 2 t — ap 2 z i t 3 + abpz 2 t + abp 2 z i t 3 — abpz 2 

pz 2 t 

+ {bt ~ b ~ t) ( 1 ~ apz2) cl>(t,p\ -pz 2 t 2 ) -(bt-b- t)^t,p 2 , -p 3 z 2 t 2 ) . (4.16) 
pzH 

To obtain an explicit expression for Q w one first expands the g-product in the function (f> with 
the help of the (/-binomial theorem [33] to obtain 



oo oo 



</>(t, q, xt 2 ) = 1 + t m Q n{n ~ 1] 

m=0 ra=l 

where the (/-binomial coefficient is defined as 



m — n — 1 
n — 1 



(4.17) 



N 



(<?;<?) 



N 



< M < N 



(q;q)M(q;q)N-M 

otherwise. 



(4.18) 



After some algebraic manipulations this gives us an explicit and surprisingly elegant expres- 
sion for Q w as 



Q w (z,a,b,p) = ^2(-pz 2 ) m p 2m ^x 

m=0 



((1 - b) 


w — m 


+ b 


w — m + 1 




m 


V 2 


m 



Cl-a)(l-6) 



w — m + 1 


- (l-a)b 


w — m 




m — 1 


m — 1 




p 2 





(4.19) 



and hence via (|4.10p an explicit expression for the generation function for loops Lwi 
L w (z,a,b,p) = 



E (-pz 2 ) m p 2rn ^ m ~ 1 ^ (1-6) ["-"»] + b[ 

m=0 V 



w— m+l"| 
m 



E (-p^pMm-i) ((1 - 6) [»--] p2 + 6[--+ 1 ] p2 - (1 - a)(l - 6) r-^] p2 - (1 - a)6[»-_?] p2 ) 

(4.20) 



m=0 



Note that the sums involved only have finitely many non-zero terms, as the (/-binomial coefficients 
are zero when m > [w/2 + lj . 

For a = b = 1 we obtain the particularly simple identity 



L w (z, 1, l,p) = 1 



E (- P z 2 ) V m(m_1) ["] 



m=0 



m=0 



(4.21) 
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5 Analysis of the free energy 



5.1 Finite width 

When the width is finite (w < oo) the system is effectively one-dimensional and so there cannot 
be any phase transitions at finite temperatures. Mathematically, one can see this in the following 
way. Regardless of whether p = 1 or p < 1 the generating function L w (z,a,b,p) is a ratio of 
polynomials P w (z,a,b,p) and Q w (z,a,b,p). As such the only singularities in the generating 
function occur at zeros of Q w (z, a, b,p). Given that these are orthogonal polynomials satisfying 
the same three-term recurrence, albeit with different initial conditions, the zeros of P w (z, a, b,p) 
are distinct to those of Q w (z,a,b,p). The simple zero in Q w leads to a simple pole in L w . The 
location of the zero z c (w; a, b,p) is an analytic function of a, b and p. 

5.2 Half-plane limit 

We now use that 



lim 

w— >co 



w — m 
m 



_ p 2 (p 2 ,P 2 )m 



to give 



lim Q w (z,a,b,p) = V(-pz 2 )"V 

w— >oo * » 

m=0 



2m(m— 1) 



(1-a)- 1 



(5.1) 



(5.2) 



and 



{p 2 ,P 2 )m (p 2 ,P 2 )rn-l 

= H(z,p)-(a-l)pz 2 H(p 2 z,p) 
= aH(z,p) + (1 - a)H(pz,p) 



lim P w (z,a, b,p) = lim Q w (z, 0, b,p) = H(pz,p) . (5.3) 

to— >CO lO-^OO 



One can therefore demonstrate explicitly that for p < 1 

lim L w (z,a,b,p) = L hp (a,p) = —— H< f. Z,P \ u( r 

w^oo aH(z,p) + (1 - a)H(pz,p) 



(5.4) 



The series H(z,p) converges absolutely for any z provided p < 1, and so has no singularities 
as a function of z. Therefore, the only singularity of L hp (a,p) occurs as a result of zeros of 
H(z,p) — (a — l)pz 2 H(p 2 z,p) provided they don't cancel with zeros of the numerator. As such, 
consider that the denominator can be written as H(pz,p)—apz 2 H(p 2 z,p) and that the numerator 
is H(pz,p). For small z we have that H(z,p) is positive since H(0,p) = 1. Consider the smallest 
positive zero of the H(pz,p): let us call it z p . Now since p < 1 we have H(p 2 z p ,p) > so 
that the first zero of H(pz,p) — apz 2 H(p 2 z,p) must occur at some value z c < z p : that is, no 
cancellation occurs. 

Necessarily the closely singularity to the origin of L hp (a,p) is then an analytic function of a. 
Hence there is no phase transition as a function of a. 
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The situation for p = 1 is different and well-known. It has been shown [10] previously that 

hm L w (z, a, b, 1) = L hp {a, 1) = . (5.5) 

1 + VI - 4z 2 - 2az 2 

Considered as power series in z, the denominator (and numerator) of the half-plane generating 
function only converges for z < 1/2. For small o the algebraic singularity at z = 1/2 is the 
closest to the origin while for a > 2 the pole arising from the zero in the denominator is closer. 
This leads to a phase transition on varying a. 



5.3 Infinite slit limit 

Let us consider the limit for p < 1 of 

lim z c (w; a, b,p) . (5-6) 

As we have discussed above z c (w; a,b,p) arises from a zero of the polynomial Q w (z,a,b,p) 
while Zc P (a,p) arises from a zero of H(z,p) — (a — l)pz 2 H(p 2 z,p) which is analytic for all z. Now 
we have that the limit of Q w (z, a,b,p) is the half plane denominator H(z,p) — (a — l)pz 2 H(p 2 z,p) 
for all z. Hence, since the zero at Zc P (a,p) does not cancel with a zero in the numerator as shown 
above we can deduce 

lim z c (w; a,b,p) = z^ p (a,p) . (5-7) 

w— >oo 

We also know that z c (w; a, b,p) is an analytic function of a and b for any fixed w and p and that 
Zc P (a,p) is also analytic in a. 

Importantly, the argument fails when p = 1 because the limit of the denominators for finite 
widths have a factor, depending on b, that does cancel with one in the numerator. 



6 Comparison of p < 1 with p = 1 cases 



To summarise for p < 1 we have just argued that the free energy of our model is the same in 
the infinite slit and the half plane scenarios: 



K inf - Slit (a,b,p) = lim K(w;a,b,p) = K hp (a,p) 



(6.1) 



for all a, b and that K hp (a,p) is an analytic function of a. On the other hand we know |10| that 
for p = 1 



K 



inf—slit 



(a, b, 1) = lim k(w; a, b, 1) = < 



' log(2) if a, b < 2 

log (-t=t) if a > 2 and a > b (6.2) 



log 



otherwise. 



while that 



n hp (a, 1) 



log(2) a < 2 



(6.3) 
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Hence for b > max(a, 2) 

K inf - Slit {a,b,l) / K hp (a,l) (6.4) 

and so the two-wall infinite slit scenario is physically different to the one wall half-plane for 
p = l. 

From a physical point of view the p = 1 case is more complicated than the p < 1 cases in that 
there are phase transitions in both the half-plane and infinite slit scenarios and, crucially, these 
are not all coincident. For p < 1 the infinite slit and the half-plane are the same. This implies 
that for p < 1 (negative pressure) in the infinite slit the top wall plays no part physically in 
the behaviour of the system. For p = 1, while the walls are a macroscopic distance apart in the 
infinite slit, the polymer can still feel both walls, recalling that the polymer is also macroscopic 
in length. 

We have solved a simple model of a vesicle confined in a slit. Mathematically the solution is 
interesting as it provides a new orthogonal polynomial series associated with Dyck paths. Phys- 
ically it demonstrates that an osmotic pressure term can remove the effect of confinement and 
the effect of a contact potential with a far wall no matter how strong the potential. Interesting 
further work would be to investigate the scaling around the limit w — > oo and p — > 1 . 
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